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Abstract 

We continue to study the spherically symmetric vacuum solutions 
of the tensor-four-scalars theory which is a modification of general rel- 
ativity. Our aim is to construct vacuum solutions with asymptotically 
constant circular velocity curves which can be used to represent the 
dark halo outside the baryonic matter. These solutions are obtained 
by means of a new non-standard gauge in the field equations. The 
metric is expressed in terms of the rotation curve V(r) in closed form. 
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1 Introduction 



As pointed out in the first part of this series [1] the tensor- four-scalars theory 
is the classical theory corresponding to the massless limit of massive gravity. 
On the theoretical side, taking quantum gauge theories for spin-1 and spin- 
2 gauge fields as the basic framework to describe nature [2-5], this theory 
has the same right for being considered as fundamental as Einstein's theory. 
In this respect our theory differs from the many modifications of general 
relativity discussed in the literature. On the experimental and observational 
side the theory is very interesting because, due to the four massless scalar 
fields, it has additional degrees of freedom which tentatively might explain 
the dark matter phenomenology without introducing hypothetical particles. 

There is a second important aspect of the dark matter problem: the right 
choice of the gauge. Since gravity is described by a spin-2 gauge theory 
different gauges are possible. In the case of static spherically-symmetric 
problems described by means of spherical coordinates t, r, <f> the standard 
gauge defined by 522 = 9m = —r 2 is widely used. The corresponding vacuum 
solution in general relativity (disregarding rotation) is the Schwarzschild 
metric. It is unique due to Birkhoff's theorem. This choice of gauge is 
appropriate for small bodies up to stars and clusters, but obviously not for 
galaxies. The reason is that the corresponding circular velocity (rotation 
curve) must decrease as V 2 ~ 1/r. Until now no such decrease has been 
observed in galaxies. Therefore, it is highly desirable to investigate other 
gauges which lead to solutions with different asymptotic behavior. This is 
one purpose of this paper, the other being the tensor-four-scalars theory. 

The paper is organized as follows. In the next section we start from a 
discussion of the measurement of the rotation curves in galaxies. This will 
lead us to the introduction of a new non-standard gauge. In section 3 we 
study spherically symmetric vacuum solutions in this gauge. We reduce the 
field equations to a single second order ordinary differential equation in the 
radial coordinate r. In section 4 we integrate this equation and express the 
metric functions by the circular velocity V(r) in closed form. In the last 
section we fit the tail of the rotation curve of M33 given by Corbelli [8] to a 
simple vacuum solution. 
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2 Choosing a non-standard gauge 



Dark matter has been introduced in astronomy to explain the rotation curves 
of galaxies. Therefore, we must discuss how the circular velocities and the 
radial distances are measured. We use spherical coordinates = (t, r, (f) 
throughout which define our laboratory system on earth. The radial distance 
r is measured in light years or kpc. The distances of nearby galaxies can 
be determined by measuring the apparent luminosities I of some standard 
candle (as the classical Cepheids) with a priori known absolute luminosity 
L. Then the so-called luminosity distance is given by 




In the large we assume Euclidean geometry, the space is only locally curved 
in the vicinity of galaxies. On the other hand we only consider nearby 
galaxies so that the cosmic expansion can be neglected. The radial distances 
r within the galaxy measured on earth as angles a in arcmin are converted 
to kpc by 

r = adL dreo- 

The rotation velocity of stars or neutral hydrogen gas is measured by 
means of the Doppler shift of spectral lines 

!^l = l-V ml + 0(V 2 ). (2.2) 

Here v is the frequency of the light source at rest and f Q bs the frequency 
actually measured by the observer on earth, V Te \ is the component of the 
relative velocity V along the direction from observer to the light source. 
Following Weinberg ([6], chapter 3/2) we introduce the freely falling coor- 
dinate system £ a of the moving star. In this system the star is at rest and 
the observer on earth moves with the 4-velocity 

u a = ^ = (u°,&), u° = (1 - V 2 Y l l\ u = u°V, (2.3) 
ar 

because in the locally inertial coordinates special relativity holds. Here r is 
the proper time 

dr 2 = Va ^ a d^, (2.4) 
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Va/3 = diag(l, —1,-1,-1) is the Minkowski tensor and the speed of light is 
put = 1. The Doppler shift is given by the special relativistic result ([6], 
eq.(2.2.2)) 

~ - TTk7' (2 - 5) 

which yields (2.2) in leading order. 

In addition Weinberg introduces a laboratory coordinate system x^ which 
in our case is attached to the observers telescope. To simplify the following 
discussion we assume that the astronomer on earth has corrected his mea- 
surements for the motion of the earth with respect to the center of the galaxy, 
so that we can choose the center of the galaxy as origin of the laboratory 
coordinate system. Now the star moves on a geodesic ([6], eq.(3.2.3)) 

d?x^ dx a dx 13 

^r + C— — = (2.6) 

where F^a are the Christoffel symbols of the metric g^ u . The latter is defined 
by ([6], eq.(3.2.7)) 

Since the local inertial coordinates £ a can be changed by arbitrary Lorentz 
transformations, we can choose a Lorentz boost such that the two coordinate 
systems are at rest with respect to each other. Then we have 

£ = 0, and ^ = (2.8, 
for j = 1,2,3. This leads to 

/0£°\2 A d? dC . , 

500 = teJ ' 9 * = -±<WdxS (2 ' 9) 

and goj = 0. Now it follows 

n d£,° d£° dx° , dx° , n , n . 

" " IF ~ WTF - {2M > 

In the following we assume that the non-diagonal elements gjk vanish. 
Then from the invariant 

a B , CK2 -.2 d?*_ d^dx^dx" 



u a ri a pu p = (u ) 2 - u 2 = 1 = -rr—rlap 



dx^ '" p dx" dr dr 
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dx^ dx v /cix u \2 /dx J \2 

= 9 ^^ =m y-d7) +9 "VdT) 

we find 

T?2 _ 1 „ fdx^2_ * 9jj (dxi/dT) 2 

v ~ ( u oy9n{ dT )- ^ gm{dx o/ dT) 2 ^ iL > 

where (2.10) is used. 

We want to specialize this for circular motion r = const. Note that in 
this case there is no first order Doppler effect for an observer in the center 
of the galaxy because V Te \ = 0. But an observer on earth has V ve \ / if he 
sees a spiral galaxy under some angle of inclination i ^ 0, for example. Now 
the geodesic equation (2.6) for a = 1 reads 

$ + ^)Vh(£)V^)Vri(^o. (, 12 ) 

Taking the circular orbit r = const in the equatorial plane •& = ir/2, then 
(2.12) gets simplified to 

itr-W 

and from (2.11) we finally obtain the important relation 

V 2 = (2.14) 

ffoo r 33 

In the following we consider static spherically symmetric metrics of the 
form 

ds 2 = g^dx^dx" = e a dt 2 - e b dr 2 - rV(cM 2 + sin 2 tf# 2 ), (2.15) 

where a(r), b(r), c(r) are functions of r only. These metric functions are so- 
lutions of the gravitational field equations (Einstein's equation or the tensor 
4-scalars theory). It is well known that the field equations do not determine 
the g^ v uniquely there is freedom to choose a gauge. In the standard gauge 
one puts c = 0. It is argued that this can be achieved by redefining the 
radial variable as 

r' = re c/2 . (2.16) 

However, this new r' is no longer the observable radial distance as defined 
above. Furthermore, c(r) as solution of the field equations is not a pri- 
ori known. Consequently, there is no clear definition of r'. Even if c(r) 



5 



were known, it would be different for different galaxies, so that r' would be 
different; a very unpleasant situation. 

The strongest argument against the standard gauge c = comes from 
observations. The vacuum solution with c = is unique (Birkhoff 's theorem) 
and from (2.14) one finds V 2 ~ 1/r for large r. This does not fit the 
measured rotation curves for galaxies which show V — > Vfl a t = const. The 
standard way out is to abandon the vacuum equations and assume some 
hypothetical dark matter. As long as this dark matter is not convincingly 
recorded one should study the other possibility of retaining c(r) / 0. Then 
the vacuum solution is no longer unique (see Sect. 2 below). To fix it uniquely 
we take the expression (2.14) f or the circular velocity V(r) as our non- 
standard gauge condition. In other words, (i) every circular velocity V(r) 
defines a particular gauge, (ii) since the circular velocity is an observable 
it is gauge invariant in contrast to the metric, therefore, our gauge fixing 
has a direct physical meaning. However, that means V[r) must be given, 
so that the theory seems to have less predictive power. But what seems to 
be a weakness is a strength: V(r) cannot be predicted on the basis of the 
vacuum equations alone, the dynamics of the normal matter must necessarily 
be taken into account. Indeed a uniform asymptotic velocity profile seems 
not to exist in reality. There is the baryonic Tully-Fisher relation [7-8] 
which relates the asymptotic velocity Vfl at to the fourth root of the total 
baryonic mass M 1 / 4 . But beside this the velocity curves at large r have 
quite different behavior for different types of galaxies. Our gauge condition 
(2.14) precisely allows for this. The non-standard theory is much richer than 
standard general relativity where asymptotic flatness is usually assumed. 

It has been argued that the function c(r) specifies the relation between 
the radius r and the circumference of circles, and therefore c(r) should be 
zero. This geometric interpretation of c(r) can be used for small bodies 
where both r and the circumference can be measured by comparison with 
a measuring rod. But in a big galaxy one cannot do so. Then the (gauge- 
dependent) metric is only accessible by observing gravitational effects as the 
circular velocity (2.14). An additional geometric interpretation has no phys- 
ical justification. We share Weinberg's non-geometric views ([6], sect. 6/9): 

"The geometric interpretation of the theory of gravitation has dwindled 
to a mere analogy, which lingers in our language in terms like 'metric', 'affine 
connection', and 'curvature', but is not otherwise very useful. The important 
thing is to be able to make predictions about images on the astronomers' 
photographic plates, frequencies of spectral lines, and so on,..." 

The physical correctness or failure of our theory can only be decided by 
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working out its further consequences (and by observing dark matter particles 
of course). Concerning this we refer to part III of this series which is now 
available in the arXiv 1205.4309. 

3 Spherically symmetric vacuum solutions 

According to [1] the vacuum field equations of the tensor-four-scalars theory 
are the modified Einstein equations 

~ \a^R = ^^{-d^v n d u v n + \g^g aP d a v n d p v n } (3.1) 

and the wave equation in the metric g a ^ 

d a (V=jjg afi df } v n ) = 0. (3.2) 

Here the index n is raised and lowered with the Minkowski metric, but all 
other Greek indices with g^ v as in ordinary general relativity. G is Newton's 
constant. 

To study static spherically symmetric solutions we choose spherical co- 
ordinates 

x ° = t, x 1 = r, x 2 = x 3 = 4> (3.3) 

and set the speed of light c = 1. We assume the metric to be of the following 
non-standard form (2.12) 

5oo = e a , gn = -e b , 522 = -r 2 e c , g 33 = sin 2 $g 2 2 (3.4) 

and zero otherwise. Since we look for static solutions we neglect the time 
derivatives do in the wave operator (3.2): 

pd/2-c 

0a(V=O/ a %)) = -smW r (e d / 2 ~ b r 2 d r )-e d / 2 - c d^sm^) - —^d 2 

(3.5) 

where d = a + b + 2c. Then the wave equation for a static v n becomes 

^f e d/2-6 r 2^W = e d/2-c L 2 v n (3 6) 

or V or J 

where L 2 is the quantum mechanical angular momentum operator squared. 
Since we restrict ourselves to spherically symmetric solutions the right-hand 
side vanishes. Then only the time component v n = (vo(r), 0, 0, 0) can be 
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different from 0, otherwise we would get non-diagonal elements in the metric. 
The equation (3.6) can now be integrated once 

v'o(r) = ^ 2 e b - d/2 , (3.7) 

where A is an integration constant. The prime always means partial deriva- 
tive with respect to r. 

The Christoffel symbols and the Ricci tensor for the metric (3.4) have 
been computed in [1]. Only the diagonal elements are different from 
zero. Using (3.7) the modified Einstein equations (3.1) for the 00, 11, 22 
and also 33 elements of the Ricci tensor give the following three differential 
equations 

c" = + h'c> + V - 3c') + ^ - 1) + (3.8) 

= ia'c' + V + c') + ^ + l(l - e b ~ c )+^je 2b ~ d (3.9) 



2 r 4 r 2 V / 2r 4 

a" + c " = -(&' - a' - 2c') - -(a' 2 - a'b' + a'c' - b'c' + c' 2 ) + - A e 2h ~\ (3.10) 
r 2 r 4 

where 

a = 16ttGA 2 (3.11) 

is another form of the constant of integration for v' (r). In the following it 
is convenient to subtract (3.8) from (3.10) 



'' ! = d -T~ t («' -b' + c') + -(c' - a') + 1(1 - ( > 
4 2 r r 2 



+ (3.12) 

and to use this equation instead of (3.10). 

It is not hard to eliminate the function b(r) from the equations. We solve 
eq.(3.9) for b(r) 

b = c + 2 log r - log(l - ^ e_a_C ) + 

. /a'c' a' + c' c' 2 1\ , 010 , 

This enables us to eliminate 6 in the remaining two equations (3.8) and 
(3.12). As we have already mentioned in [1] there is a degeneracy in these 
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equations. Since this is an important property we give some details of the 
elimination process. Differentiating (3.13) we get 

, 2 / 2r 2 \ / , , 2\ 

b' = c' + - + (l - — 2 -) (a' + c' + -) + 

r \ 2r z — a exp(— a — c) ' \ r' 

+^ (a"(2r 2 c' + 4r) + c"(2r 2 a' + 4r + 2r 2 c') - 4(a' + c + ^)) , (3.14) 
with 

N = 2r 2 a'c + 4r(a' + c) + r 2 c' 2 + 4. (3.15) 
We substitute this into (3.12) and solve for c": 



2 lr' <■'- 1 x jV 



( ' + + + 1 
a' + c' + 2/r \ r ra' a' r 2 a' ' 2r 2 — a exp(— a — c) r 2 

(3.16) 

This enables us to simplify (3.14) considerably: 

b' = 2^ + a' + 2c , + -. (3.17) 
a' r 

Next we turn to (3.8) and substitute (3.17) for b' . This yields 

" a " ( * i ^. 2 i 2rVc / + 4r(q / + c / )+rV 2 +4 
c =— (c+-)+^H T1 : : • (3.18) 

This is identical with (3.16) which shows the degeneracy in the system of 
vacuum field equations. Summing up from (3.18) we find a class of solutions 
with one free function a(r) or c(r). This reflects the freedom of gauge. 

4 Integration of the vacuum equations 

It is convenient to simplify the equations by introducing the new metric 
function 

f(r) = c(r) + 2 log- (4.1) 

r c 

for c(r). The quantity r c has been included for dimensional reason. Then 
eq.(3.18) assumes the simple form 
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and the circular velocity squared (2.11) simply becomes 



V 2 (r) = j = u. (4.3) 

There should be no confusion with the 4-velocity u a in (2.3) because this 
quantity only appear in Sect. 2. We also introduce the new variable 

a + C I A A\ 



so that 



. 2u 1 + rw 

a = 

r 1 + u 

, 2rvJ - u 

c = — . (4.5) 

r 1+u y ' 



Then (3.18) becomes 



u' 2r 1 + i\ n 

— + 77^ ? — 7^ — \~i 1 + rw ) = °- 

u 2r z — aexp(— 2w) 1 + u 

After multiplication with (l + u)/(l + 2u) this is integrable in the form 

Hence we get the integral 

u 2 {r 2 e 2w -a/2) _ u 2 f o a+c 



l + 2u l + 2u 



r< "f)=/5. (4-7) 



P is a new constant of integration. 

Next we have to incorporate the circular velocity (4.3) which leads to 



2 

r 

On the other hand from (4.7) we get 

a+c 1 / „ 1 + 2u . a 



a' = u(c' + -). (4.8) 



e 



Taking the logarithm and differentiating we obtain a' + c'. Then using (4.8) 
we separately calculate 

2?/' 

a' = -— — (4.10) 

1 + 2u + 7n 2 y ' 
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with 



u(l + 2u + 7ii 2 ) r 

7 



20" 



(4.12) 



Now everything is reduced to simple quadratures. The solution depending 
upon the existence of real roots of 1 + 2u + ju 2 , we set, for 7 < 1, 



71 = 1 + x/1 - 7, 72 = 1 - V 1 - 7- 
Then that we arrive at 



(4.13) 



v 1— 7 



a = < 



7 

it+1 
2 



+ 



1+72^ 



I+71 it 

og K a 



+ logK a 



arctan 



Ay— 1m 
l+u 



+ log K a 



for 


7 


< 1 


for 


7 


= 1 


for 


7 


> 1 



(4.14) 



with K a an integration constant. Here we have chosen the indefinite integral 
in such a way that it is continuous over < u < 00 and — 00 < 7 < 00. In 
particular, near low velocities (remember that u(r) is the circular velocity 
squared (4.3)) 



a = —2u 



l-u + 



^-^u 2 + (^-2)u 3 + 0(u 4 )} + logK a V 7 , 



whereas near 7=1 
2u 



l + u 



1+ ^(T^) 2(1_7,+ KTT^) 4(1 " 7,2+0((1 " 7 ' 3) ] +logii ' ' 



and near 7 = 0, which is the special case of ordinary general relativity, 
a = -log(l + 2u) + 



- \ log(l + 2u)j 7 + 0(7 2 ) + log K a . 
1 + 2u 2 J 



Function c(r) obtains immediately from a(r) using (4.12) 



c = —a + log 



1 1 + 2u + ju' 

r 2y2 



+ log K a K c 



(4.15) 



with K c another integration constant. Substituting these results into (4.9) 
we get the following relation between the integration constants: 



K a if c = /3sgn(l + 2u + 7u 2 ). 



(4.16) 
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K„ 



1+72" 
1+71 u 



K a e "+i 



arctan ■ 



k K a e ^ 
In the same way we find from (4.15) 



l+u 



Finally b(r) also follows from a by (3.13), (4.7) and (4.15) 

b = -a + log {^f + log K a K c 
Taking the exponential we get the 00-element of the metric 

1A/1=7 

for 7 < 1 

for 7=1 

for 7 > 1 

for 7 < 1 
for 7 = 1 
for 7 > 1 

for 7 < 1 
for 7 = 1 
for 7 > 1 



(4.17) 



(4.18) 



K, 



|l+2u+7« 2 | 



l+7i u 
1+72 u 

, 1+2 + w\ ^ 

K |l+2«+7» 2 | — — 
-<V. Z7T771 e 



7-1 



arctan ■ 



(4.19) 



l+u 



and from (4.17) 



e = < 



K c 

K r 



l+7i m 

1+72 M 
2 2u 

1 



l/v^ 



(4.20) 



e v^FT arctan 



/7— 1m 
l+u 



As a check, we can specialize these results to the classical Schwarzschild 
case which must be contained, of course. To this end, we first specialize our 
result to the pure tensor theory of general relativity, i.e. to the case a = 0. 
Then 7 = by (4.12), and 71 = 2 and 72 = by (4.13), so that (4.18-20) 
becomes 

K a h /u'\2 (l + 2u) 2 

e a = — V, e b = K c (-) l + 2« , e c = K C K J . 4.21 

Next we assume the Schwarzschild metric in standard form, as given eg. by 
[6] (8.2.12). By (2.14) the circular velocity squared is then equal to 

r2 r s 



V 



u = 



(4.22) 



2(r - r s ) 

with r s = 2GM the Schwarzschild radius. Inserting (4.22) in (4.21) yields 

r e< = ^. (4.23) 
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After identifying 

K a = 1, K c = ^ (4.24) 
we recover the Schwarzschild metric and the standard gauge c = : 




(4.25) 



But clearly it would be dogmatic to forbid non-standard gauges. Moreover, 
in order to match the outer vacuum solution to an inner solution with normal 
matter, general relativity (4.21) might not be flexible enough. To be on the 
safe side we therefore continue to investigate the tensor-four-scalars theory. 



5 Fitting the tail of the rotation curve of M33 

Since we have neglected baryonic matter until now, we can compare the 
vacuum solution only with the tail of a real rotation curve where the con- 
tribution of the visible matter is small compared to the dark. Probably the 
best object for such a comparison is M33. This is a low luminosity spiral 
in the Local Group, it is dark matter dominated and near to us. There are 
fairly good data available which have been analyzed by Corbelli [8]. We 
are indebted to Edvige Corbelli for communicating the numerical values of 
Fig. 5 of his paper. 

As we have shown here the rotation curve V(r) is not restricted by the 
vacuum solution. However our experience with power series expansions in 
[1] suggests that for large r a series in powers of r _1 should be a good 
representation of 

V\r) = Vl t + ^ + ^ + ... = u(r) (5.1) 

According to the analysis of Corbelli the visible and dark contributions to 
V 2 (r) are equal around 7 kpc. Therefore, we have taken the last 8 points 
with r > 7 kpc of the data and fitted them with weights given by a to four 
terms of (5.1). As the figure shows we obtain a very good fit with a reduced 
X 2 = 0.1868. This well passes a \ 2 - test with 99% confidence level and 
8 — 4 = 4 degrees of freedom. Of course this is not a test of the theory. A 
real test can only be carried out if additional data beside the rotation curve 
are available, for example lensing data. We have discussed such a test in [5]. 
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Finally we note a peculiarity of the metric corresponding to the asymp- 
totic behavior (5.1). From (4.16) and (4.18) we find the asymptotics of the 
metric functions 

- gn =e b = ^ + 0(r- 5 ) (5.2) 

-522 = rV = L 2 + O^ 1 ) (5.3) 
and similarly for 533. This gives for the determinant of the spatial part 

-D = -511(722533 = e b+2c r 4 sm 2 $ = ^iM sin tf + 0(r" 5 ). (5.4) 

Since the volume measure is proportional to \J — D it follows that \/—D ~ 
r~ 2 is integrable at infinity, so that the metric gives a finite volume to 
the 3-space (r, cf>). That means the metric can only be locally true in a 
galaxy and its near neighborhood. In the large V 2 (r) must deviate from 
(5.1). Regarding M33 Corbelli has found that the "dark halo" extends to 
a distance comparable with the separation distance between M33 and its 
bright companion M31, the Andromeda galaxy. The model of one isolated 
galaxy in infinite space seems to be an unrealistic idealization. Therefore, 
it is no harm that we have no theory for this at present. 

Acknowledgment. We thank the referee for his criticism. It has en- 
abled us to improve the manuscript considerably. 
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Figure caption: Figure: M33 rotation curve (points) with a fit to 
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